Geometric Manipulation of Ensembles of Atoms on Atom Chip for Quantum 

Computation 
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We propose a feasible scheme to achieve quantum computation based on geometric manipulation of 
ensembles of atoms, and analyze it for neutral rubidium atoms magnetically trapped in planoconcave 
microcavities on an atom chip. The geometric operations are accomplished by optical excitation of 
a single atom into a Rydberg state in a constant electric field. Strong dipole-dipole interactions and 
incident lasers drive the dark state of the atom ensembles to undergo cyclic evolutions that realize 
a universal set of quantum gates. Such geometric manipulation turns out naturally to protect the 
qubits from the errors induced by non-uniform laser illumination as well as cavity loss. The gate 
performance and decoherence processes are assessed by numerical simulation. 

PACS numbers: 03.65. Vf, 03.67.Lx, 42.50.Pq, 32.80.Ee 



I. INTRODUCTION 

Quantum information processing (QIP) holds out the 
possibility to run algorithms and protocols superior to 
those of its classical counterpart [HQ. In recent decades, 
numcrious candidates for physical implementation of 
quantum information processors have been proposed Q . 
Because of long coherence times and exceptional con- 
trolability, quantum optical and atomic systems such as 
trapped ions d, Q , neutral atoms [1, Q , and cavity QED 
[1,0] , have taken a leading role in implementing quantum 
logic. The photon is a remarkably robust candidate for 
a qubit, and easy to transport long distances; however, 
as we know, they interact only weakly with each other, 
which makes the realization of quantum gates based on 
photons difficult. For instance, in linear optical quan- 
tum computation(LOQC) (lOj . extra measurements are 
required. On the other hand, ensembles of trapped atoms 
or molecules may serve as convenient and robust quan- 
tum memories for photons. They can act as an interface 
with flying qubits, and store and retrieve single photons, 
for example by electromagnetically induced transparency 
(EIT) [ll|, [l2| , which allows to keep coherence of quan- 
tum states up to several hundreds of millisecond |l3[ . 

One might therefore hope to use controlled interac- 
tions in the ensemble of atoms, once the flying qubits are 
stored, to realize universal quantum logic gates in a de- 
terministic and scalable way. Several promising schemes 
for such operations based on combining EIT with Ryd- 
berg atoms was proposed in Ref. fl4l . fl5| , exciting atoms 
with lasers to high-lying Rydberg states and exploiting 
the strong long-range dipole-dipole interaction between 
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Rydberg states. Recent theoretical proposals have ex- 
tended towards a single step high-fidelity entanglement 
of a mcsoscopic number of atoms [l|| and quantum sim- 
ulation both coherent dynamics and dissipative evolu- 
tion process for many body system (l7j . Remarkably, 
the building blocks behind all such setup-up have been 
demonstrated experimentally by several groups [H, [l!| . 
In Ref. [2(| , an effective strong dipole-dipole blockade ef- 
fect is achieved by coupling to microwave coplanar waveg- 
uide resonators to realize quantum logic on ensembles of 
molecules. Here, we put forward an alternative, hope- 
fully scalable approach to realize universal quantum gates 
based solely on laser-controlled geometric manipulations 
of neutral atomic clouds trapped on the surface of atom 
chips. The use of geometric gates can naturally eliminate 
certain kinds of errors. 

The paper is organized as follows. In section II we 
introduce the basic model of atoms and cavity system. 
Then, we show how to construct the Hamiltonian of the 
system and the method of adiabatically manipulating the 
system evolution to implement three quantum gates ge- 
ometrically: the phase gate, the y-rotation gate, and the 
controlled phase gate. These constructions are given in 
sections IIIA, IIIB, and IIIC , respectively. Numerical 
analysis of their performance is also given in these three 
sections. These three gates forms a well-known univer- 
sal set of quantum logical gates. Finally, in section IV, 
we discuss how to perform qubit measurements on this 
system, and conclude in section V. 



II. MODEL 

Atom chips supply a good platform on which precise 
control and manipulation of the neutral atoms can be 
performed [2l|, [22|. We will show that an atom chip 
with integrated Fabry-Pcrot (FP) microcavity can real- 
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FIG. 1: Schematic of trapped ensembles of atoms in magnetic 
atom traps inside plano-concave optical microcavities (hori- 
zontal). Each atomic cloud is a qubit, and forms a processing 
cell. Different qubits couple to each other mediated by the 
photonic modes of another Fabry-Perot cavity (vertical). 



ize a universal set of all-optical control universal quantum 
gates for atomic ensemble qubits. 

Consider several plano-concave Fabry-Perot (FP) mi- 
crocavity resonators [23| integrated on the atom chip, 
which are parallel to each other. Another large FP cavity 
is integrated so that its mode is perpendicular to all the 
plano-concave cavities. The diagram is shown in Fig. [T] 
The plano-concave resonator consists of an isotropically 
etched dip in a silicon surface, and the cleaved tip of a 
single-mode fibre 2J] serves as the input-output channel. 
Atomic ensembles are placed in the region of highest field 
strength of the cavity modes, leading to higher values of 
the coupling g, and a Q value over 10 6 . The perpen- 
dicular FP cavity modes serve as a data bus to couple 
between different qubits. 

The atomic clouds must be confined in traps inside 
the cavities. We can introduce magnetic fields produced 
by current-carrying wires on the surface and coupled to 
the magnetic dipole of the atoms. The small scale of 
the structure produces strong magnetic field gradients, 
which make tight traps for magnetic atoms. Atoms in 
a weak-field-seeking state will be attracted and held in 
this region. The hyperfine states \F = 2,mf = 1) and 
\F = 1, mp = —1) of 5S 1; /2 of 87 Rb are ideal weak-field- 
seeking states that can be trapped by a static magnetic 
potential [2l|. In addition, these two states have op- 
posite Lande factors, so that they experience identical 
magnetic potentials. Thus, the decoherence of an arbi- 
trary superposition of these two states due to current 
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intensity fluctuations is strongly inhibited 
oscillations between these states have been observed with 
decoherence time as long as r c = 2.8 ± 1.6s [2(|. So, it is 
natural to consider these two states of 87 Rb for quantum 
information processing. 

The state of the 87 Rb atom cloud can be initialized 
by inputting a single photon qubit to the plano-concave 
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FIG. 2: Storage of single photon qubit. Level scheme of a 
single 87 Rb atom is shown. The input qubit is encoded as 
jO) if there is no photon in the input mode and |1) if there 
is a photon in the input mode. Here, we use the hyperfine 
levels of manifold 5Si/2: clock states \g) = \F = 1, rriF = 
-1) and \s) = \F = 2,m F = 1); |a_) = \F = l,m F = 1) 
and \a+) = \F = 2,m,F = —1) are used as ancillary states. 
And a state of manifold 5Pi/2 serves as an intermediate state 
|e) (\F = l,rriF = 0)). The storage process mapping single 
photon signals into collective excited atomic states of \s) can 
be accomplished by adiabatically turning on the classical laser 
field Q(t) to a value much larger than the coupling constant 
g when the photon arrives. The reverse process can be used 
to read out the photon. 



cavity through the fibre, as shown in Fig. [2j The state 
|0) p h corresponds to no photons in the input channel, and 
the state |l) p h corresponds to one photon. An arbitrary 
qubit state can be represented as cv|0) p h + /3|l) p h- In 
general, the state of n input channels is an entangled 
state in a Hilbert space of dimension 2". 

An ensemble of N identical multi-state atoms is 
trapped in each cell. Using well-developed techniques, 
all atoms can be initially prepared and trapped in a spe- 
cific sub-level (hyperfine level \gi), i = 1, 2, 3, N) of 
the internal atomic ground space of manifold hj. Rel- 
evant states of each atom include the other clock state 
|sj), two ancillary states |a_^}, |a+j), and an intermediate 
state \ei). (As mentioned earlier, |sj) has a long coher- 
ence time.) We also assume the atomic density is not too 
high, so that interaction between atoms can be safely ne- 
glected. Atoms are manipulated by illuminating the en- 
tire ensemble, to excite all atoms with equal probability, 
so only symmetric collective states are involved in this 
process [141 ] . In the case of loading quantum information 
into the ensembles, we consider only three states: |ej), 
| Si) and \gi), where |e*) and \gi) are coupled to the cavity 
mode of the plano-concave cavity, while |e,) and \si) are 
coupled by a classical field. Initially, all the atoms are in 
their ground states: 



|g) = |g 
We define operators 
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The storage state with excitation n is 



|g w -",s") 



(N -n)\N r > 
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For the case when n = 1, this is 
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|s) = -f= 92, 9s-9n) + \gi,s 2 , gs-gN) 



H 1" \9i,92,93--sn))- 

The commutation relation of iS and S' is 

For a sufficiently small number of excitations (n <C 2V), 
these two operators can be approximated by bosonic an- 
nihilation and creation operators of quasiparticles corre- 
sponding to \s). Operators £ and £\ and A± and A± 
have the same properties. 

As shown in Ref. and Ref. [28[, when the pho- 
ton arrives in the cavity, one can adiabatically turn on 
the classical field coupling the states |e) and |s), until 
the Rabi frequency f2 is much larger than the cavity- 
atom coupling constant g, so that the state of photon 
is stored in the ensemble of atoms in the cavity in the 
form |ri)ph — > |s"), where |n) p h represents the n-photon 
Fock state. This process is reversible: one can retrieve 
the photon from the ensemble of atoms. So, an arbi- 
trary input state a|0) p h + /3|l) p h can be mapped into the 
state a|g) +/3|s) with high fidelity H3- Then, We denote 
|0)l = |g) and = |s) for the logical representation of 
a qubit in each cloud of atoms. 



III. A UNIVERSAL SET OF GEOMETRIC 
GATES 

We now discuss how to geometrically implement quan- 
tum computation (holonomic quantum computation (29| ) 
in our system. The first step is to construct a Hamilto- 
nian that has an eigenspace with eigenvalue to avoid dy- 
namic phase during the cyclic evolution. This eigenspace 
can either be nondegenerate, which would introduce sim- 
ple Abelian phase factors (Berry phases) [3(| or denerate, 
which could cause general non- Abelian operations [3l| . It 
is well known that single-qubit operations, together with 
a nontrivial two-bit gate, make a universal set of quantum 
logic gates for quantum computation [32j ■ By construct- 
ing an appropriate set of looped paths in the parameter 
space we can obtain an arbitrary unitary transformation 
in the computational space. 



We choose for our universal set of gates 
J?«(0i) = exp Uh\U) LL {\ 



b 2 ) = cxp i(j>2<J. 



exp 



^3|lj0fc) i L<l i 0fc 



The states |0j)i and \li)L are the computational basis 
states for qubit i, and a\ is the Pauli operator in the y 
direction for the ith qubit. It is well know that if we can 
implement these gates with arbitrary 4>\ , 4>y. , 03 , we can 
implement all unitary transformations [321 ] - The three 
operators correspond to a rotation in the z direction, a 
rotation in the y direction, and a conditional phase rota- 
tion when qubits j and k are in the state |10)l. In this 
section, we show how to holonomically realize these three 
gates. 



A. One Qubit Phase Gate 
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FIG. 3: Left: Level structure and laser coupling diagram for 
a single atom. Right: The equivalent coupling diagram for an 
ensemble of atoms in second quantized representation. 

We first consider how to perform a simple one qubit 
phase gate: — > e l ^' 1 \l)L- To simplify the model, we 
only consider the energy levels we are interested in. The 
total pulse sequence and resonant coupling diagram is 
shown in Fig. [3] The Hamiltonian of of this system, in 
the interaction picture and rotating wave approximation, 
is 



ff 1 (t)=^oi(t)(^i ) +H. c ) 

i 



(1) 



T< jk. = are transition operators from 

state \k) to \j) for atom i, and h = 1. Here, the two 
lasers are assumed to be phase matched, and their Rabi 
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frequencies are real numbers. For simplicity, assuming 
the laser beams illuminate each atom with same intensity, 
we can rewrite the Hamiltonian in the second quantized 
representation: 



basis, we have 



+ n a {t)^A+ + n* a {t)Als, 



(2) 



and set fli = fisinfl and O a = — ficosfle 1 ^. The pa- 
rameters 6 and ip are functions of time, and the absolute 
magnitude Q is a constant that must be large enough to 
satisfy the adiabatic condition. £,S, A+ are bosonic op- 
erators defined in last section corresponding to the single 
atom states |e), \s), \ a+), respectively. 
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FIG. 4: (a) The evolution of the relative phase between state 
and state |0)l for three different values of f2, to real- 
ize a 7r/I6 gate. Note that after the cyclic evolution, we 
get an additional pure geometric phase, which fits the the- 
oretical value quite well for f2 larger than 200MHz. (b) The 
population change of state during the process. We can 
see that only if f2 is large enough to satisfy the adiabatic 
condition can we get high accuracy of the gate. Two short 
laser pulses are determined by 8(t) — exp((t — 1) 2 /0.I5) and 
ip(t) = exp((t - 1.59I7) 2 /0.f5) (the unit of t is /is). 



We can regard the system as three coupled harmonic 
oscillators. We are only interested in singly excited 
states; the Hamiltonian discussed above is closed in 
that basis, which for simplicity we represent as {|e) = 



1 



|s) = \N - 1,0,1,0), 



\N - l,l,0,0) se 
|a + ) = ]7V — 1, 0, 0, l) g esa}- The coupling diagram of the 
ensemble of atoms is equivalent to a simple three-state 
coupling diagram for a single atom, as shown in Fig. [3] 

Representing the Hamiltonian in matrix form in this 



Hx(t) 



flsmO 
n sin e 
-ficosfle-^ 








(3) 



The zero-energy eigenspace of the Hamiltonian is non- 
degenerate. The dark state is 



\D{t)) = cos0|l) L + sin6>e^ ¥ '|a + ). 



(4) 



Note that the Hamiltonian is decoupled from the state 
|0> L = |A,0,0,0) grsQ . 

Assume the initial state is oo)) = \D(— oo)) 
1 ; (which means 8(—oo) = at the very beginning). 
Also, assume that initially ip = 0. Now, turn on two 
laser beams, using the standard formula to calculate the 
geometric phase, where the parameters make a cyclic evo- 
lution with the starting and ending point 6 = 0, ip = 0: 



0i = i I dR(D(t)\V R \D(t)}, 

where R(i) = (0(t),ip(t)). We have 

(D(t)|V R |D(t)) =-ism 2 (e)<p 

and 

sin 2 (6')d^. 



(5) 



According to Green's theorem, this geometric phase is 
exactly the enclosed solid angle: 



shi9d0dip = dfi. 

OS JJdS 

A pulse sequence consisting of two stimulated Raman 
adiabatic passage (STIRAP) pulses with relatively large 
width (to guarantee that the process is adiabatic) can be 
used for this purpose. 

A numerical simulation of resonance coupling to real- 
ize a 7r/16 gate (^i = ir/8) is shown in Fig. |U Since 
in the laboratory the Rabi frequency 57 can be as high 
as 200 MHz, the adiabatic process could be finished in 
several microseconds, which is much shorter than the de- 
coherence time of the hyperfine ground states. We can 
exchange the population of and |a+) back and forth 
without losing any information, and we can get a geo- 
metric phase with high accuracy if fl is large enough. 

In the previous method of using ensembles of atoms or 
molecules to realize quantum information processing, a 
major shortcoming is that it is very hard to uniformly il- 
luminate all the atoms, thus limiting the accuracy of the 
gate operation. But in our proposal, each atoms interacts 
individually with Ox and O a , and the geometric phase is 
independent of the absolute value of f2. If the spatial 
distribution of the electric fieldsjof the two jaser beams 
match each other, say f?i(r) = Ex(r)d es = E a (r)d ea+ = 
Q a (r) for every point in space r, the effect of non-uniform 
illumination can be naturally eliminated. This is an ad- 
vantage of adiabatic control. 
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B. Y-Rotation Gate 



l4>2 &y 



Now we show how to achieve the gate U = e 
This gate is much more complicated to realize than the 
one qubit phase gate. It is difficult to couple the state 
|0)l to since directly coupling them by applying 

a laser beam to the ensemble of atoms would result in 
unwanted higher collective excitations. In Ref. [3], a 
scheme taking advantage of the dipole-dipole interaction 
of Rydberg states between atoms was used to realize this 
gate. Here we propose a geometric gate by following three 
steps: 

1. Adiabatically pump the ground state |0)l to a 
highly excited Rydberg state |r), using the dipole- 
dipole interaction to assure that only a single col- 
lective excitation is achieved. 

2. Adiabatically control the coupling between the sin- 
gle excitation states |r), |s), and |a_) to geometri- 
cally realize the gate. 

3. Reverse step 1 to transfer |r) back to |0)l. 

To carry out step 1, consider Rydberg states of a hy- 
drogen atom within a manifold of fixed principle quan- 
tum number n with degeneracy n 2 . This degeneracy can 
be removed by applying a constant electric field E along 
a certain axis (linear Stark effect), such as the z axis. 
For electric fields below the Ingris- Teller limit, the mix- 
ing of adjacent n manifolds can be ignored, and energy 
levels are split according to AE nqm = inqeasE /2 with 
parabolic and magnetic quantum numbers q and m, re- 
spectively. Here, q can take values n — 1 — |m|,n — 3 — 
\m\, — (n — 1 — |to|), e is the electron charge, and as 
the Bohr radius. These states have dipole moments of 
jjL = 3nqeaBe z /2. 

Consider two atoms k and I seperated by a distance 
R. The dipole-dipole interaction between them is 



^fp(R) 



i 



47T60 



IRP 



,(£* • R)(£i • R) 



IRI 5 



(6) 



where fl k and jli are dipole moment operators for atoms k 
and I. Suppose the electric field is sufficiently large that 
the energy splitting between two adjacent Stark states is 
much larger than the dipole-dipole interaction. 

For two atoms initialized in Stark eigenstates, the di- 
agonal terms of V^p(R) provide an energy shift, while 
the nondiagonal terms couple adjacent m manifolds with 
each other, (m, m) — > (m ± 1, m =F 1). The off-diagonal 
transition can cause decoherence, and is suppressed by 
an appropriate choice of initial Stark eigenstate Q ■ (For 
hydrogen, the state is \r)„ = \n,q = n — l,m = 0).) For 
a fixed distance R = Re 2 , the dipole-dipole interaction 
of two atoms k and / in the states \rk) ni and |n)n 2 is 

<„ 2 (R) = nM\nM\Vgmr k ) ni \n)n 2 





Stpl 



FIG. 5: The level diagram with the couplings to realize the 
Ry{4>) gate. We apply an electric field to the ensemble of 
atoms, and pick appropriate Stark eigenstates of rubidium to 
be the register states \r) = |r)7o,o, |/) = k)72,o and the inter- 
mediate Rydberg state \m) = |r)60,o as shown the Appendix. 
States \g), |a_) and \s) are the in the same manifold of the 
internal ground state. In the first step, two conjugate laser 
beams fi p i and Q P 2 are used to adiabatically pump a single ex- 
citation from |g) to |r). In the second step, three laser beams 
fio, fii, are used to adiabatically realize the R x y {<f)) gate. 
The excitation from \g) to \r) is then adiabatically reversed. 



For ni and sufficiently large, u(R) cx nfn|. For al- 
kali atoms such as rubidium, the situation is more com- 
plicated, but this kind of analysis still works. In the 
appendix, we discuss how to pick a state which has char- 
acteristics similar to \r) n for rubidium. We will make use 
of this energy shift later. 

Consider the diagram shown in Fig.[SJ For step 1, there 
should be a single excitation in state |r) after pumping 
if the initial state is \0)l, and no excitation in state |r) if 
the initial state is To achieve this goal, we introduce 
two other Rydberg states |m) and |/), and adiabatically 
transfer the amplitude of to |f) by coupling \s) and 
|/) to some intermediate state, say |e). Please note that 
this pumping process can also be realized adiabatically. 

The Hamiltonian of the cloud of atoms as described in 
the left side of Fig. [S] in the rotating wave approximation 
is: 



H p (t) = ^2(Q pl (t)afl + H.c) 

i 

+ ^(O p2 (t)a« +H.c) 



k>l 



u mm (R-ki)\m k )\mi)(m k \(mi\ 

-u^(R fci )k fc )|n)(r fc |(n|) 

-uf f CR kl )\fk)\fi){fk\{fi\ 
- u^ m (R k i)\r k )\mi)(r k \( mi \ 
-uf m (R kl )\f k )\ mi )(f k \( mi \ 



(7) 
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The two-body interaction term here is due to the energy 
shift of the states. Just as in the case of the phase gate, 
we represent the Hamiltonian in second quantized form. 
Again, assuming the laser fields are uniformly coupled to 
each atom, and using the standard second quantization 
procedure to deal with the two-body interaction, we ap- 
ply the minimum energy shift for all pairs of atoms (as 
the worst case) and obtain 

iCl pl (t)VN(M^ +M) 

+ U mm M ] M ] MM + Ur^^TUl (8) 



(*) 



u ff- 

V>fm 



where u rr , u m m, w//, u rmi Uf m and Hf r are chosen to be 
the minimum energy shifts from the dipole-dipole inter- 
action. 



N-l, 1,0,1 




|0) =|N,0,0,0) 



I N-l, 1,0,1) 




Ufl - Ufm 



N-l, 1,0,1 



FIG. 6: The equivalent coupling diagram of the ensembles of 
atoms in cases of collective excitation for the case (a) system 
is initially in \0)l and (b) system is initially in 

In the case of low excitation, in the Bogoliubov ap- 
proximation we can regard the operators 0,0^ on state 
\g) as the C-numbcr \ where TV is the total number 
of the atoms in the cloud. When the system is initially in 
\0)l (and hence |f) is not excited), the equivalent state 
coupling diagram for this pumping process is shown in 
Fig. EJa). Unlike the phase gate, we see that this Hamil- 
tonian is no longer closed in a small subspace of states 
because fl p i will continuously pump atoms of the ensem- 
ble from the ground state to the excited states. How- 
ever, because of the large energy shift, the rate to excite 
|iV-2, 1,0, l) srsm and |iV-2, 1, l,0) grsm is strongly sup- 
pressed. Physically, this means there cannot be an atom 
excited to the Rydberg state |r) while another atom be 
excited to state \i) or \m). Only a single excitation to 
the register state can be achieved (|0)l — > |r)) during 
this adiabatic process. 

If the system is initially in the state |1)l (and hence 
there is a single excitation in |/)), both |r) and \m) can- 
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FIG. 7: Adiabatic Pumping process. Note that at the end 
there is only a single excitation in the register state |r) when 
initially in |0)l, and no excitation when initially in Here, 
we set JL rm = 400MHz and 7I mm 300MHz. Since we take the 
minimum energy shift as our simulation parameter, the ac- 
tual performance should be better (i.e., the process could be 
finished in a shorter period of time). 



not be excited due to the energy shift caused by the 
dipole-dipole interaction, as shown in Fig. [Bfb) . A nu- 
merical simulation of the pumping process for the mini- 
mum energy shifts is shown in Fig. [7J 

After the pumping process, and the adiabatic transfer 
of the population of |f) back to |1)l, we can then realize 
the Ry(4>) gate on the ensemble of atoms by applying 
three controlling laser beams. The coupling diagram is 
shown in the right side of Fig. [5] The Hamiltonian of step 
2 can be obtained directly in second-quantized form: 



+ {l a (t)(£Al+&A-). 



(9) 



We can directly apply the technique of geometric trans- 
formations to realize this gate, as introduced in Ref. [33| 
for four-level systems. Here, we choose Slo = ^ sin 9 cos <p, 
Hi = £ls'm8 simp, Q a = £lcos(9. This time, the Hamil- 
tonian is closed in the basis {|e),|0)z, = |r) , 1 1) z, = 
|s),|a_)}: 



H 2 (t) = n 



o 

sin 9 cos ip 
sin 9 sin ip 
cos 9 



sin v cos ip sin u sm p cos ( 



(10) 

The eigenspace corresponding to the zero-energy eigen- 
value (dark space) is spanned by basis vectors {|Di), and 
\D 2 }}, where 

\Di) = s\np\Q) L - cos^|l) L , 



\D 2 



cosf9(cos(^|0)L + sin<p|l)i) 
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FIG. 8: (a): The pulse shape of the control lasers of Slo, ^i and Q a respectively. (b)(c)(d): The evolution of the populations 
of \s), |e), \r). We see that when fi is larger than 200MHz, the adiabatic condition is satisfied and the gate works with very 
high accuracy (close to the theoretical value) and the population of the short-lived state |e) tends to zero. 



We can use the formula for the degenerate subspace 
under the cyclic evolution of 9, (p. Suppose \4>(t)) = 
Ci{t)\Di) + C 2 {t)\D 2 ). We have the equation 



(11) 



"Ci" 






012 (*) " 










_D 21 (t) 


022 (t) . 




c 2 



where Dij(t) = (Di\Dj). Let's define 
dA 2 (t) = 



011 (t) 012 (*) 
021 (t) 022 (*) 



— cos 8tp 
cos 6 ip 



and also the unitary matrix 



U =Vexp(i <j> dA 2 (t)). 



After the cyclic evolution, 

|Di(T)} =W 11 |0 1 (O))+W 21 |0 2 (O)), 
|0 a (T)) =W 12 |0 1 (O))+W 22 | J D 2 (O)). 



(12) 



(13) 



(14) 



At the beginning, if we set 6(0) = and (p(0) = |, we 
will have |0i(O)) = |0) L and |D 2 (0) = \1) L . We can use 
the Dyson expansion to get the unitary matrix U. This 



gives us: 
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u 
u 

u 21 
u 22 



l _ il + Tl 



2 24 720 

= sin ^2, 



cos <p 2 , 



_ 4 if. 

6 120' 

b-2 h . . . = — sin (t> 2 , 

2 6 120 v ' 



(15) 



24 720 



cos <p 2 , 



where 



4>2 = J cos Odip 

is a pure geometric phase. Thus we get the unitary trans- 
formation we want. Step 3 is simply the reverse of step 
1, to return us to our original basis. 

A simulation to realize <p 2 = 7r/4 was done for a cer- 
tain pulse shape of the coupling laser beams, where the 
qubit was initially prepared in the state |0) £,, as shown 
in Fig. |SJa). We show the evolution of the populations 
of states |s), |e) and |r) in Fig. [5fb),(c), and (d), re- 
spectively. If f2 is large enough (> 200 MHz, which is 
practical in current experiments), we could implement 
the gate operation in Ifis with extremely high accuracy. 
Thus, the total procedure takes less than 3/xs to com- 
plete, while the lifetime of |r) is approximately equal to 
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Atom in qubit 1 










a 




g+ u 






t A 







Atom in qubit 2 



FIG. 9: The schematic setup to realize the conditional phase 
gate, and the relevant atomic levels of the two clouds of atoms. 



the lifetime of |r) [jjj . which is estimated to be 300-400/is 
in a cryogenic environment [34| . 

We simplified the calculation and simulation by as- 
suming that all atoms are uniformly illuminated in the 
laser beam. However, just like the argument made in the 
previous section, the geometric phase is independent for 
large f2, if the spatial distribution of the electric fields 
of the three laser beams match each other at each point 
in space. Thus, we can in principal eliminate the error 
introduced by nonuniform illumination. 



C. Two Qubit Controlled Phase Gate 

Now, we are ready to construct the controlled phase 
gate. Recently, an adiabatic SWAP operation between 
states of two clouds of atoms in a cavity QED system 
was experimentally achieved (35j . Here, we extend this 
method to geometrically realize the controlled phase gate, 
which together with the earlier one-qubit gates gives a 
universal set of quantum gates. The basic idea is to 
couple two qubits (that is, clouds of atoms) by virtu- 
ally emitting and absorbing a cavity photon; and after 
an adiabatic evolution, the second qubit obtains an ex- 
tra relative phase when the first qubit is in \0)l- The 
coupling diagram of the scheme is shown in Fig. [9] 

Suppose the length of the cavity is about 200/xm, which 
means, it could contain about 10 qubits. Consider two 
different clouds of atoms in the cavity. For a single atom 
in each cloud, pick a manifold of states of p. When an 
electric field is applied to these two cloud of atoms, the 
manifold splits into |r+) and |r_). For atom cloud 1, we 
set a constant electrical field E z \ and a constant magnetic 
field B z i in the z direction. The Stark Effect splits the 
manifold p, and the Zceman effect changes the energy of 
\g) and \s). 



Assume the energy difference between \g) and \s) to be 
(5. Set the appropriate non-zero value of E z \ and B z \ so 
that |r + ) is resonantly coupled to \g) and \s) by a cavity 
mode a and a control laser Oi respectively. For atom 
cloud 2, E Z 2 is set so that \r + ) is resonantly coupled to 
\g) and |s) by the same cavity mode a and control laser 
0,2, and |r_) is resonantly coupled to |s) and |a_) by the 
cavity mode and a control laser Q, 3 . Note that for Rb 87 , 
the hyperfine structure energy split of the ground state 
manifold is A/27r=6.835GHz, and the energy difference 
of |r + ) and |r_) should set to be equal to A for atom 
cloud 2 for the purpose of resonant coupling, while the 
frequencies of fl\ , 0,2 and 3 should be cu cav — 5, ui cav — A 
and uJca V +A, where ui cav is the frequency of cavity mode. 

Since the cavity mode is inhomogcncously distributed, 
the thermal motion (even at extremely low temperature) 
of atoms, causes the coupling rate between atoms and 
cavity photons to vary greatly from one to another (by 
roughly a factor of 2) and thus makes the system diffi- 
cult to control accurately. To overcome these difficulties, 
instead of directly applying the laser beam to the clouds 
of atoms, we use the idea of external laser driving control 
[36| as shown in schematic setup in Fig. 9. Three classi- 
cal laser fields E\ (t) , E2 (t) and E 3 (t) are incident on one 
mirror of the cavity to drive the transition \r + ) — > \s) for 
cloud 1, and \r + ) — > \a + ) and |r_) — > |a_) for cloud 2 
through another cavity mode a! . We assume for simplic- 
ity that a and a' have the same spatial mode structure 
x( r ) with the same frequency uj C av, so they can be differ- 
ent only in polarization. 

The driving fields E\ (t) , E2 (t) and £"3 (t) are resonant 
with frequencies u; cav — S, uj cav — A and uj cav + A. So 
they are far-off-resonant to the cavity mode with large 
detuning. E\{t), E 2 (t) and E 3 (t) control the time evolu- 
tion of the Rabi frequencies fii, O2 and f^. To see this, 
consider the input-output equation for the cavity mode 
a' [13: 

= -i[a(t),H sys ] - -a'(t) + v^(t), (16) 

where a' ia (t) is the field operator for the input driving 
pulse coupling to the mode a', with 

(a[ n (t)) = E^t) + E 2 (t) + E 3 (t), 



and 



L(t),4(t')]=6(t-t>). 



Since a' is driven by strong classical pulses, we can treat 
the mode as classical and thus assume that the interac- 
tion between a' and the atoms will not change the state 
of the cavity mode a' . Eq. (fT6| therefore can be modified 
to be 



da'jt) 
dt 



iuj cav a'(t) - 2 a 'i t ) 
sfH{Ex{t) + E 2 {t) + E 3 (t)). 



(17) 
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First, we suppose there is only one classical pulse E\(t) 



-S)t 



, where £\{t) is the slowly 



varying amplitude of E\ (t). Then, we can write the mean 
value of a' as 

(a'(t)) = a(i)e-^*- 5)t , 
which has the solution 

ri 



a(t) 



£i(t) exp(— iS — n/2)(t — r)dr 
id + k/2 



(18) 



when the characteristic time T\ of £i(i) satisfies 5T\ 3> 1. 

Assume that e±(t) gradually increases from zero with 
£i(0) = 0. Thus we can conclude a(t) cx e±(t). By choos- 
ing an appropriate phase of £i(i), we set a(t) to be real. 
If we input E2(t) or E 3 (t), the solution will have a similar 
form. Now, consider the case when three classical pulses 
are incident on the cavity mirror. Eq. (|16p is a linear 
equation, so, when the three input pulse have different 
resonant frequencies, the solution should be a superposi- 
tion of the solutions for three single pulses. Thus, (a'(i)) 
could be represented as 



(a'(t)> = ai(i)e 



i ( u 



-8)t 



-a 3 (t)e~ 



+ a 2 (t)e~ 

+A)t 



i ( u 



-A)* 



This gives Rabi frequencies fij(r, t) = di<a!j(t)x(r) for i = 
1, 2, 3, where di is a coefficient mainly determined by the 
dipole moment for the corresponding transition. 

The frequencies of the three components of a! differ sig- 
nificantly, so we can regard them as three separate pulses. 
In addition, the energy structure of atoms in cloud 1 is 
different from those in cloud 2, so ill will not interact 
with the atoms in cloud 2. Similarly, fi 2 and O3 will 
not interact with atoms in cloud 1. All other clouds in 
the cavity are all far off resonant from both the cavity 
mode and the control lasers, and thus can be safely ex- 
cluded by our gate operation. Now, we can represent the 
Hamiltonian for the situation of of Fig. [5] in the interac- 
tion picture: 

H 3 (t) = £ n x ( rjl , i)^) + £ g + (r n )aV£> 
h h 

h 

32 

+ H.c, 

where a and cr are annihilation and creation operators of 
cavity mode a, and g+{v) and <7-(r) are coupling rates of 
cavity mode a to the transition \g) |r+) in cloud 1 and 
\s) —> |r_) in cloud 2. We can represent g+(r)(g- (r)) = 



(19) 



g + x(r)(5_x(r)), where 17+ and g_ are constants. So we 
can rewrite the Hamiltonian as: 



31 

+ E X(r n )(d 2 a 2 (t)a^l + + d 3 a 3 (t)a^± 
31 

+ ff + aV^+S-aV«)+H.c., 



(20) 



We will see soon that the dark state of such a Hamil- 
tonian should be independent of x( r )j an d the adiabatic 
evolution of the system is determined only by the mini- 
mum value of x( r j) f° r au h sa Y Xm- So, the Hamiltonian 
can effectively be represented as: 



Hf(t) = ^(niCty&o + 5+a t o -«0) 
+^(o 2 (^) + + o 3 (t)^L) 



(21) 



+ V(ff+ata^ + ff_aVW) + H.c 



t„C»ah 



ja 

where fi»(i) = di<Xi{t)xm f° r * = 1,2,3 and g+{g~) = 
g+Xm(g-Xm)- As in the previous cases, we transform 
the Hamiltonian to second quantized representation. For 
two clouds of atoms in the Bogoliubov approximation, 
we get the new Hamiltonian 

Hf(t) a n 1 (t)(^) 1 (TZ+) 1 + n 2 (t)(AlMn+h 

+ a 3 (t)(Al)2(n.) 2 + g + VNa^(n + ) 1 

+ g+ VNcS(K+)2 + g-a t (5 t ) 2 (7t_) 2 + H.c, 

(22) 

where S, 1Z-, and A+ are the bosonic annihi- 

lation operators of the single atom states \s), |r_), \r + ), 
|a„), and |o+), respectively. We set 

fi 3 (i) =fi 3 cos 6e~ llp3 . 

We also define g\ = g + vJV, g 2 = g+ \/N and g 3 = g- for 
convenience. 

As in the previous sections, for each cloud of atoms we 
have a set of basis states: 



|0>l = 


\N, 0,0,0, 0,0) gr+r _ a+a 


- s 1 




|1>l = 


|JV-l,0,0,0,0,l) sr+r _ 


a_|_a- 


-S) 


r+> = 


|A-l,l,0,0,0,0) gr+r _ 


a_|_a_ 


-Sj 


r-) = 


|^-l,0,l,0,0,0) 9r+r _ 


a_|_a_ 


-Si 


a+> = 


|JV-l,0,0,l,0,0) 9r+r _ 


a_i_a_ 


-Sj 


a_> = 


|JV-l,0,0,0,l,0) 9r+r _ 


a+a_ 


. s ■ 
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The Hamiltonian is closed in a subspace of Hilbcrt space 
that can be divided into the direct sum of two closed 
subspaces: 



(23) 



Here, H 1 (E)'H 2 _ is spanned by basis vectors { 1 100) , |r+00), 
|001), |0r+0), |0a+0)}, and U 1 ®U 2 _ is spanned by basis 
vectors {|110>, |r+10), |011), |0r_0), |0a_0)}. The first 
two degrees of freedom represent cloud 1 and cloud 2, and 
the last degree of freedom of the state is the Fock state 
of the cavity photons. Since the Hamiltonian doesn't 
couple these two subspaces, and has the same form in 
each subspace, we can just consider a single subspace, 
e.g., U 1 ®U\. The Hamiltonian can be represented in 
this subspace as 



H 



eff 



(*) = 




9i 


Sli sin# 






5i 

O2 cos 9( 




91 
9i 








fi 2 cos 6e l f2 






Qi sin 9 







(24) 



where we have set (pi = 0. A dark state exists for this 
system, since one eigenstate has eigenvalue 0: 



\D(t)) = 



Tm cos0 



nffe cos2 6 + sin2 1 



= |100) 



cos^ f sin 



Ittt COS 2 



n?(t 



cos 9 sin ( 



cos^ f sin 



Ittt C 



cos^ f sm 



|0a + 0) 



1001) 



(25) 



Just like the case of the phase gate, if the system is ini- 
tially prepared in the dark state, after a cycle of adiabatic 
evolution a pure geometric phase can be obtained. In this 
case, if the initial state is |10)l, we get a phase shift 

<t> 3 = i I dR(D(t)\Vn\D(t)) 



dip 2 



92 



9l 



(26) 



Note that this phase is not affected by the spatially in- 
homogeneous distribution of the cavity mode. 

We also get a geometric phase for the subspace V} ® 
H 2 .. However, we can set choose our path to set this 
phase to independently of <p^. So, if we initially pre- 
pared the state |10)i we get a geometric phase of 03 after 
the gate manipulations, and if the initial state is 




FIG. 10: The evolution of the phase 4> 3 of |10)l. The pulse 
shapes of 9(t) and ip(t) are also given. 



the geometric phase is 0. The states |00)l and |01)l are 
not affected by this Hamiltonian, and hence also acquire 
no phase. The net effect is a two-qubit quantum gate, 

|00) L ^ |00>£, 
|01>£ -> |01) L , 
|10) L ^e lfe |10) L , 
|11>l^|H)l, 

which is exactly the controlled phase gate we would like 
to implement. 

As we can see in Eq. (|2"5]) , the two clouds of atoms inter- 
act with each other by virtually absorbing and emitting a 
photon in cavity mode a, so leakage of cavity photons is 
the most important source of decohcrence in our scheme. 
Analytically, the average photon number in the cavity 
during the process is 



cos 2 9 sin 2 9 



at 



92 

am 



sin 9 + cos 2 9 sin 9 



(27) 



when the system is initially prepared in |10)l, and 



cos 2 9 sin 2 1 



T-ph- 
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-P- 



(28) 



when the system is initially prepared in |11)l, respec- 
tively. If the initial state is |00)l or |01)l there are no 
photons in the cavity. 

Let's consider the concrete example of a conditional 
phase gate with ^3 = n/16. We set the parameters to 
be g+ = 20MHz, g 3 = g_ = 10MHz, and the number of 
atoms in each cloud to be N ~ 10 3 . Then we have g\ = 
660MHz, fii = 40MHz, fl 2 = 50MHz and Q 3 = 300MHz. 
We choose 



9(t) 



■exp[-(t- 1.7) 2 /0.5], 
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Va (f) = ~exp[-(t-2.115)70.5], 

(where time is expressed in /is). A numerical simulation 
(not including cavity loss) is shown in Fig. [10] where the 
initial state is \\00) L + \\01) L + ±|10) L + J|11) L . The 
whole process is very fast and can be finished in 4/zs, 
with extremely high accuracy approaching the theoretical 
value calculated in Eq. ([2"6)) . 

Fig. [TT] shows the average number of cavity photons 
n p h+ and n p h- during the process as a function of time 
and ill (not including cavity loss). The photon number 
is less than 0.001 in general which means that, the prob- 
ability of a photon loss in the cavity is bounded above 
by 0.001 for the parameters given previously. We see 
that when the state of system is initially in the 
average number of photons in the cavity is larger than in 
1 10) i, |00)i and |01)l. (Indeed, for the last two states 
the photon number is strictly zero.) 




0.5 1 2 4 8 



k(MHz) 

FIG. 12: Fidelity of the controlled phase gate with cj>3 — it/ 16 
for different values of k and g\. 

Suppose the cavity loss rate is k. If that is the domi- 
nant decohcrence channel, the fidelity of our conditional 



gate can be given by 0: 

F = min \/(^\m(TT/16)C(\^)(^\)U(Tr/16)\^}, (29) 

W 

where the superoperator C is determined by the master 
equation: 

g = -i[Hf, p] + 2n[apa^ a, p}} (30) 

in the Markov approximation. As stated previously, the 
initial state |11)l emits more photons than the other 
three basis states during the gate operation. Hence, we 
can take \ip) = as the worst case. 

From Eq. (|26[) , we can see that is approximately in- 
dependent of (?i and when they are large enough, so we 
can compare the fidelity for different value of g\ for the 
same operation. Fig. [T2l shows the gate fidelity for differ- 
ent values of k and g\ in our gate operation. The simu- 
lation is realized with CH — h library "Quantum State Dif- 
fusion" [HI using a fourth-order Runge-Kutta integrator 
[39l ] and a pseudo-random number generator to solve the 
master equation using quantum trajectories. (For each k, 
we computed only 15 trajectories, so this is not a highly 
precise calculation.) It makes sense that g\ = 1000MHz 
has better performance than that of 171 = 660MHz, since 
the large value of g\ decreases the number of cavity pho- 
tons and hence gives a higher fidelity. The fidelity can 
exceed 0.999 when k — 500kHz for gi = 660MHz, and 
when k = 1MHz for g Y = 1000MHz. For a FP cavity 
larger than 100/im, values of k under 1MHz are quite 
achievable. 



IV. MEASUREMENT 

In fault-tolerant quantum computation, one must do 
syndrome measurement periodically @, l40l l4l| , and the 
measurement results should be correct with high proba- 
bility. It is also necessary to read out results at the end of 
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a computation. Accurate measurement is therefore nec- 
essary. In this section, we propose a protocol to measure 
an atomic ensemble qubit in the computational basis. 

Collecting the fluorescence from an atom is the most 
natural way to realize a measurement. For a single atom 
(or ion), many fluorescence photons are needed to make 
the measurement reliable, so a "cycling" transition from 
the computational state (|0)l or to an unstable ex- 
cited state is generally used. In our case, however, there 
is an added complication. The computational states in- 
clude superpositions of a single excitation over all the 
atoms of the cloud. The distance between atoms in the 
cloud can be several times larger than the wavelength of 
a photon emitted from the cloud, and this might make 
it possible (in principle) to distinguish which atom emit- 
ted the photon. This in turn could cause decoherence 
by collapsing the symmetric superposition, and take the 
system state outside the computational space. So we 
need a method to perform a reliable measurement while 
avoiding this problem. 




Is) l a ") Is) l a -) 

Stpl Stp2 

FIG. 13: The schematic setup to realize a projective measure- 
ment on a single qubit. Here, state \d) is a temporary state 
with a short lifetime. The transition between \d) and \s) is 
not allowed. The state \d) can be chosen for example, to be 
\F — 2, uif = —2) of the manifold 5P3/2. 

As shown in Fig. 13, the procedure can be realized in 
three steps: 

1. Adiabatically transport the population of |1)l to 
|f). (The procedure is similar to that used in real- 
izing the gate R y (</>).) 

2. Dynamically pump the atom cloud from the ground 
state to state |r) with a ir pulse D.q, and subse- 
quently apply another ir pulse Qi between the state 
|r) and \d). Cycle this transition hundreds of times 
and collect the spontaneously emitted photons from 
the atom cloud. 

3. Reverse step 1 to transfer the population of |f ) back 
to to prepare for next operation. 

If the qubit is initially in |0)l, step 1 has no effect. The 
state of the cloud will transfer to |r) and then to |d) after 



two 7T pulses; the atom cloud will quickly decay to the 
ground state and emit a photon. By repeating this cycle 
many times, the probability of detecting the photons can 
be made quite high. 

If the qubit is initially in step 1 will produce a 

single excitation in |/), and by the dipole-dipolc block- 
ade effect the transition from the ground state to |r) will 
be blocked. Therefore we will not observe any emitted 
photons. Thus we observe fluorescence only if the qubit 
is initially in the state \0)l, and this procedure is a pro- 
jective measurement. Note that non-uniform coupling is 
not a problem here since at the end of the procedure the 
state will be either |0)l or 

The time for a single ir laser pulse (and the accompa- 
nying state transition) should be roughly 10-20ns when 
Rabi frequency of the pumping laser is hundreds of MHz. 
Assuming a life time of about 10ns for \d), each cycle of 
measurement can be finished in 20-30ns. After the cy- 
cling transitions, we must transfer the population of |f) 
back to |1)l, so we need the total measurement time to 
be less than the lifetime of |/) (estimated to be about 
400/is). This gives an upper limit on the total number 
of cycles we can make during the process. If we want to 
limit the error to a reasonable level, we can perform at 
most a few hundred transition cycles. This should be ad- 
equate if the photodetector has sufficiently high efficiency 
and large enough solid angle of detection. 



V. CONCLUSIONS 

In this paper we have proposed a scheme to achieve 
quantum computation using geometric manipulation of 
ensembles of atoms in a cavity QED system. Adiabatic 
optical control can be used to obtain a geometric phase 
gate and controlled phase gate. Combining optical ex- 
citation with a dipole-dipole blockade between Rydberg 
states allows us to geometrically realize the R y (4>2) gate. 
Thus a universal set of quantum gates can be realized 
geometrically. 

We analyzed this scheme for ensembles of neutral ru- 
bidium atoms, magnetically trapped in planoconcave mi- 
crocavities on an atom chip. Numerical simulations show 
that a single qubit gate can be performed in several mi- 
crosecond with very low probability of gate error. For 
the controlled phase gate, the operation is done by vir- 
tually emitting and absorbing a photon from the cavity 
mode, and can be completed in 4/is for a controlled 7r/32 
gate. An advantage of geometric manipulation is that by 
adiabatic parameter control we can avoid certain kinds 
of errors, especially those caused by inhomogeneous dis- 
tribution of the laser beam and cavity modes. The values 
of the fields can depend on atom position, but their ratio 
can be fixed and controllable. The basic idea is to find 
an appropriate adiabatic process, so that the relevant dy- 
namics are either independent of, or depend only on the 
ratio of, the two coupling rates. 

The ensemble of atoms effectively enlarges the coupling 
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rate g by y/N, which greatly suppresses the likelihood of 
cavity photons and increases the fidelity of the opera- 
tion. We analyzed the scheme for N = 10 3 atoms in each 
cloud. The parameters that we have assumed in our nu- 
merical simulations are all achievable by current exper- 
iments. This, together with the possibility of coupling 
stationary qubits (for computation) with flying qubits 
(for communication) makes this scheme look particularly 
promising for near-term quantum protocols. 
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symbols (42j : 

c(n,l,m,q) =(n,l,m\n,m,q) = (_i)C 1 -»+™+«)/a+' 

/ 71 — 1 n— 1 ^ \ 

X v2i + 1 I mXq rr3-q _ ) • 
V 2 2 m J 

(32) 

For example, if n = 2, m = 0, q = 1, we have 
c(2, 0,0,1) = \/2/2 and c(2, 1,0,1) = -\/2/2. So, 
|2,0, 1) = (V2/2)|2,0, 0) - (V2/2)|2,l,0). We will fre- 
quently use these parabolic states instead of spherical 
states in the following analysis. 

The dipole-dipole interaction Vd; p is proportional to 
r i • r 2 , where 

ri = ri( — H e y H e z ) ® I, 

n ri n 



APPENDIX 

In this appendix, we discuss how to choose the appro- 
priate Stark eigenstatcs of alkali atoms for our purpose. 
Since alkali atoms have spectra similar to hydrogen, we 
first analyze the effect of the dipole-dipole interaction of 
Rydbcrg stark eigenstates for hydrogen. After that, we 
allow for the difference between alkali atoms and hydro- 
gen, and discuss a method of calculating their energy 
structure. Using this method, we will choose appropri- 
ate Rydberg Stark eigenstates of rubidium as an exam- 
ple. (Please note that the fine structure of the p and d 
states of rubidium have observable effects on the spec- 
trum. However, for simplicity, we do not take this into 
account.) The method of calculation is described in de- 
tail in Ref. (42[. Note that all quantities below are in 
atomic units for simplicity. 

We first consider Stark states of a single hydrogen 
atom. The magnetic quantum number m is a good quan- 
tum number. From perturbation theory, the first order 
approximation of Stark eignstates are parabolic states 
\n,ni,n,2,m) [421 ] with energies 

1 3 

En,m,n 3 ,m = ~~^2 + TJ^™ 1 ~ n z) n - 

Here, n is the principal quantum number, E is the elec- 
tric field in the z direction, and ni, n 2 are non- negative 
integers satisfying the equality n = ni+n 2 + \m\ + l. For 
m = 0, allowed values of 771 — 772 are n— 1, n— 3, — n + 1 
and for m = 1, they are n — 2,n — 4, n + 2. In 
the following discussion, we use the quantum number 
q = ni — 712 instead of ni and ri2 for simplicity. 

We expand the Stark parabolic states \n,m,q) in the 
spherical basis \n,l,m) as [42j |: 



| n, m, q) = \n,l, m) (n,l,m\n, m, q). 



(31) 



r 2 = I <S> r 2 { — e x H e y H e z ), 

7-2 7-2 7 2 

are coordinate operators for atoms 1 and atom 2. Re- 
placing x, y, z with spherical coordinates r, 9, <j>, 



x = r sm a cos < 
y = r sin 9 sin c, 
z = r cos 9. 



(33) 



gives us 



Vdip oc ?i ■ r 2 

2/i , Vi 




x 2 



e z ) • r 2 ( — e x 



I An e x + ie y x 



1,-1 



' 47T e x + ie y 2 



'2 

~V2 



2/2 

— ( 

'2 



m 

7*2 



x Y 1 +Y 1 e V 

1 1.1 ' 1 1.0 e z / ) 



^ v2 



V2 



(34) 



where V^m ^ s a spherical harmonic function for atom 
i. Generally speaking, for two atoms in the given ini- 
tial Stark eigenstate, the diagonal terms of the dipole- 
dipole interaction give an energy shift, while its non- 
diagonal terms couple adjacent m manifolds with each 
other: (m, m) to (m ± 1, m =p 1). The Stark states that 
are most useful for our scheme are those that maximize 
the energy shift while suppressing the transition between 
different m manifolds (which might introduce decoher- 
ence channels). So, it is sufficient for us to know the 
value of the matrix elements 



(n, m\(n, m, q\V<u p \n' , m + l)\n' , m — 1), 
(n, m\(n, m, q\Vdi p \n', m — l)\n', m + 1), 



(35) 



1=0 



and 



The coefficients can be written in terms of Wigner 3J 



(77, m, q\(n, m, q\ Vdip I n', m, q)\n\ m, q). (36) 
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For hydrogen, we will only consider matrix elements 
with the same n for simplicity. Even so, it is too com- 
plicated (and unnecessary) to find the general analyti- 
cal form of these elements. Instead, we numerically do 
the integration to determine which states fulfill our re- 
quirements. By symmetry, the transition strengths of 
(to, to) — > (to + 1, m — 1) and (to, to) — > (to — 1, to + 1) 
are the same, so we just give the first value. The calcula- 
tion results of some matrix elements are shown in Table. U 
as an example. 

We have made several simplifications in this table. 
First, we didn't calculate the transition strengths be- 
tween different n manifolds, because they are several 
times smaller than their counterparts in the same n man- 
ifold. Second, we didn't calculate the case when two 
atoms arc initially prepared in different Stark eigenstates, 
especially in two different manifolds, which actually we 
have proposed to realize both the gate e^ 2<Ty and qubit 
measurement. Nevertheless, this table gives enough in- 
formation about the characteristics of the state we are 
looking for. First, the transition strength must be much 
smaller than the energy shift inside the same n manifold. 
Second, those transitions between initial and final states 
with a difference in parabolic number q larger than one 
must be greatly suppressed. If we prepare an atom in 
the outermost state |n, m = 0,q = n — 1), we obtain the 
largest energy shift with the smallest transition strength 
to the (m + 1,to — 1) state compared with other Stark 
eigenstates in the same manifold. 

In our proposed gate, since the Rydberg states of 
two atoms may not be in the same manifold, a natu- 
ral solution to fulfill our requirements is to choose state 
\ty) = \n, m — 0, q = n— 1) \n' , to = 0, q = n' — 1) for man- 
ifolds n and n' . Note that the energy shift term in the 
Hamiltonian is a direct product of operators on two dif- 
ferent atoms, so the analysis of two atoms in same man- 
ifold can be directly applied to show that the maximum 
energy shift is obtained if the two atoms are prepared in 
the state |\E'). This gives an energy shift of 

(*|Vdi P |*> ocn(n-lK(n'-l). 

Next, we consider the case of non-hydrogen alkali 
atoms. Physically, the main difference between alkali 
atoms and hydrogen atoms is that the former have a fi- 
nite sized ionic core that results in avoided crossings in 
the Stark spectrum. So, the Hamiltonian can be written 
as 



basis. For diagonal terms of the Hamiltonian, we have: 

1 3 

(n,m,q\H\n,m,q) = —^ + -qnE 
+ (n,m,q\V d (r)\n,m,q) + 0(E 2 ), 

and we represent the non-diagonal terms using hydrogen 
spherical states \n, l,m): 



(n, to, q\H\n', to, q') = (n, to, q\n, I, to) 

1,1' ( 39 ) 

x In, I, m\Vd(r)\n , I , to) (n , I , m\n , to, q'). 

Denote the spherically symmetric eigenstates of the 
alkali atom as \n, I, m) a i. For large n, we have 



1 



2(n-<5,) 2 



=ai{n, I, m\H\n, I, m) a i 

( V 2 1 
w(n,/,m| I — - — - + Vd(r) ) \n,l,m) 



^2 + (n,l,m\V d {r)\n,l,m). 



(40) 



where 5i is the quantum defect for angular momentum I 
of the alkali atom [43| , which is different for each element. 
For rubidium, if we neglect neglect the fine structure ef- 
fect, 5 « 3.1, 8i w 2.6, S 2 « 1.3, and 5 3 ~ 0.02. For 
I > 3, 6i ~ 0. On the other hand, 



1 



So we have: 



'2(n-Si) 2 ~ 2n 2 n 3 ' 



(n, I, m\Vd(r)\n, /, to) 



(41) 



Per Ref. [42j , this expression may be generalized to 

Si 



(n,l,m\Vd(r)\n' ,l,m) = - 



(42) 



Here, we take advantage of the spherical symmetry of 
Vd{r), so the matrix element vanishes when the I value 
on the two sides of the above equation are not the same. 

Observe that Eq. (|3l)l) can be represented as matrix 
multiplication. We define matrices 



H 



2 



+ V d (r) + Ez. 



(37) 



Here, Vd(r) is the difference of the potential function be- 
tween a hydrogen and an alkali atom, due to the finite- 
sized ionic core. We treat V d (r) as spherically symmetric 
and only nonzero near the nucleus. In the case where 
quantum defects are relatively small (e.g., n is large), we 
can use the hydrogenic parabolic states as our working 



[C nm ] q i =(n, to, q\n, I, to) = c(n, I, to, q), 
[V nm ]i q =(n, I, m\n, m, q) = c(n, I, to, q), 

[V nn '] lv ={n,l,m\Vd{r)\n',l',m), 

[S nn ' m ] qq . =(n,m,q\H\n',m,q'). 

Then, we have 

gnn hi Q-nm.'jynn. -yn m 



(43) 



(44) 
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Matrix element 


|m, g) ->■ 


\m',q') 






n =5 


n = 10 


n = 15 


n = 20 


n = 25 


|0,n 


-1) 


0,n - 


1}^ 


1, n 


-2) 


— 1, n — 


2} 


112.5 


1012.5 


3543.75 


8550 


16875 


|0,n 


-1) 


0,n - 


l>-> 


0,n 


-1) 


0,71- 


-1) 




900 


18225 


99225 


324900 


810000 


10, 71 


-1) 


0,71- 


l>-> 


1, n 


-4) 


-1, 


71 — 


4) 

















|0,ti 


-3) 


0,n- 


3}^ 


1, n 


"2> 


-1, 


71 — 


4) 137.78 


1350 


4829.32 


11769 


23362.4 


0,71 


-3) 


0,n - 


3)-» 


1, n 


-4) 


-1, 


71 — 


2) 137.78 


1350 


4829.32 


11769 


23362.4 


0,71 


-3) 


0,n - 


3}^ 


1, n 


-4) 


-1, 


71 — 


4) 168.75 


1800 


6581.25 


16200 


32343.7 


0,71 


-3) 


0,n - 


3)^ 


1, 71 


-2) 


-1, 


71 — 


2) 


112.5 


1012.5 


3543.75 


8550 


16875 


0,71 


-3> 


0,71- 


3}^ 


0,ti 


-3) 


0,71 - 


-3) 




225.01 


11025 


72900 


260100 


680625 


|1, n 


-2) 


1, n — 


2}^ 


2,71 


-3) 


0,71 - 


-3) 




137.78 


1350 


4829.32 


11769 


23362.4 


1, 71 


-2) 


1, n — 


2}^ 


1, 71 


-2) 


1, 71 - 


-2) 




506.25 


14400 


85556 


291600 


743906 


|1, n 


-4) 


1, n — 


4) -> 


2,71 


-3} 


0,71 - 


-5) 




137.78 


1350 


4829.32 


11769 


23362.4 


1, 71 


-4} 


1, n — 


4}^ 


2,71 


-5) 


0,71 - 


-3) 




168.75 


2062.16 


7744.14 


19281.9 


38742.1 


1, n 


-4) 


1, 71 — 


4}^ 


2,71 


-5) 


0,71 - 


-5) 




168.75 


2062.16 


7744.14 


19281.9 


38742.1 


1, 71 


-4) 


1, 71 — 


4}^ 


2,71 


-3) 


0,71 - 


-3) 




137.78 


1350 


4829.32 


11769 


23362.4 


1, 71 


-4) 


1, 71 — 


4}^ 


1, 71 


-4) 


1, 71 - 


-4) 




56.25 


8100 


61256.3 


230400 


620156 



TABLE I: The matrix element of operator Vdi P . 




FIG. 14: The energy level of rubidium around manifold n = 15, m = 0. The circled states is our candidate states which have 
large energy shift. 



For the purpose of illustration, we show a simple example. 
We write the matrices 



C 20 = 



V2 V? 
72 V2 



v 30 = 



V3 V3 V3 

J_ n L 

72 U V2 



1 

75 



i 



V 



23 







-JSL. 

6V6 

-6^6 



and multiply them to get 



S 230 = 



Si 



So 
36 
So 

36 T 12^/6 



Si 



So 
36 
So 
•Mi 



Sp I Si 

36 12V6 
Sp _ Si 

36 12^/e 



Matrix S nn m can be treated as a submatrix (or block) 
of the Hamiltonian matrix. Since the Hilbert space is 
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infinite-dimensional, we need to truncate the Hamilto- 
nian matrix. For example, if we need to know eigenval- 
ues and eigenstates of manifold n of rubidium, we need 
to consider only those manifolds that couple to it. 

We now calculate the energy structure of the n = 15 
manifold for the purpose of illustration. We also need to 
consider the adjacent manifolds such as n = 14, ?i = 16, 
n = 17 and n = 18 since states like 16d, 17p, 17d, 18s 
and 18p are coupled to the n = 15 manifold. Define the 
submatrix of the Hamiltonian H nn m as 

[?r m '"V = [s nn ' m W - <WV(^2 ~ \i nE )i 

where, S nm is the Kronecker delta. The Hamiltonian ma- 
trix is approximately represented as: 

r 1/14, 14, m i/14. 15, m i_/14.16.m i_/14,17,m a/14, 18, m -i 
ft ft ft ft ft 

1/1. 5, 14, m 1/15, 15, m 1_/I5,l6,m 17 15, 17, m 17 15, 18, m 
ft ft ft ft ft 

1_/I6,l4,m 17 16. 15, m 17 16.16, m 17 16, 17, m 17 16, 18, m 

ft ft ft ft ft , 

17 17, 14, m 17 17, 15, m 17 17, 16, m nj 17, 17, m 17 17, 18, m 
ft ft ft ft ft 

17l8,l4,m 1/18,15, m 17l8.l6.rn i/18,17,m a/18,18,m 
rt ft ft ft ft 

(45) 

which is an 80 x 80 Hcrmitian matrix. Similar methods 
could be applied to much higher excitation states like 
n = 80. 

By diagonalizing the Hamiltonian matrix, we can ob- 
tain the spectrum and the eigenstates. As mentioned 
earlier, the state we use should have a large component 
of \n, m = 0, q = n— 1), which means we want the single- 
atom cigenstate 

\r) n ,m = 22Cq\n,m,q) 

to have m = with a large coefficient \C q \ for q = n — 1. 

Fig. [TH shows the spectrum around the manifold n = 
15. The outermost Stark eigenstates |r)i5 j o, circled in 
the figure in the regime where the Stark eigenvalues are 



roughly linear in the electric field (E < 2.5 x 10 -7 ), are 
good candidate states, since they should have behavior 
similar to the hydrogen state \n, m = 0, q = n — 1). Nu- 
merical calculation shows that the diagonal term of the 
dipole-dipole interaction is about 8 x 10 4 in |r)i5 : o|^)i5,o 
(compared to 99225 for hydrogen in the previous table), 
which is larger than for the other states in the manifold. 

To put this in the context of our scheme for quan- 
tum computation, consider the case where there are 10 3 
atoms in a volume of (6^m) 3 . The smallest energy shift 
of a pair of atoms inside the cell is for those that are most 
distant. Suppose the two most distant atoms are in the 
state \r) i5,o\r) 15.0- The energy shift between them should 
be around 1MHz, depending on the spatial distribution 
of atoms. In our scheme, where two atoms may also be 
in different single-atom states (actually in different man- 
ifolds), we could naturally extend our original analysis, 
but we expect a similar result. 

In practice, we would want to use higher energy Ryd- 
berg states; for instance an n — 60 state \m) = |r)60,0 for 
the intermediate state and an n = 70 state \r) = |?~) 70,0 
for the register state, to get both longer lifetimes (more 
than 300-400/iS for n = 70 can be achieved in a cryo- 
genic environment [34j ]) and a stronger dipole-dipole in- 
teraction. Since higher energy Rydberg states should be- 
have more like the Rydberg states of hydrogen, we can 
pick the outermost Stark eigenstates in the linear Stark 
area for both atoms, whether or not they have the same 
principal quantum number. If the initial state of one 
atom is in the outermost state of manifold n = 60, and 
the other in either manifold n = 60 or n = 70, then 
by Vdip oc n(n — l)n'(n' — 1) the energy shift should 
be roughly 200-300MHz or 300-400MHz, respectively. 
The distance of the most closely-spaced pair of atoms 
should be about 10 4 ao- This is larger than 2R, where 
R = 4900ao is the radius of atoms in manifold n = 70. 
This spacing should maintain the validity of our dipole- 
dipole interaction model. 
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